-----------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  /users/jjacksn/documents/pathdep/conference/pdtable1.log
  log type:  text
 opened on:  15 Jun 2011, 21:18:09

. use  /users/jjacksn/documents/pathdep/conference/pdtable1

. replace y = . if[y==0]
(51 real changes made, 51 to missing)

. replace y0 = . if[y0 == 0]
(51 real changes made, 51 to missing)

. gen sample = country != 4 & country != 7 & country != 11 & country != 13 & country != 15

. replace y = y/100
(765 real changes made)

. replace y0 = y0/100
(765 real changes made)

. gen y1 = y[_n-1]
(52 missing values generated)

. gen y2 = y1[_n-1]
(53 missing values generated)

. replace y2 = y1 if[year == 1961]
(17 real changes made, 5 to missing)

. gen aussie = country == 1

. gen aust = country == 2

. gen canada = country == 3

. gen denmark = country == 4

. gen finland = country == 5

. gen france = country == 6

. gen germany = country == 7

. gen iceland = country == 8

. gen ireland = country == 9

. gen japan = country == 10

. gen nz = country == 11

. gen norway = country == 12

. gen portugal = country == 13

. gen spain = country == 14

. gen sweden = country == 15

. gen uk = country == 16

. gen us = country == 17

. gen t1 = (year - 1961)/10

. gen t11 = t1[_n-1]
(1 missing value generated)

. gen t2 = (year - 1961)/10

. gen t21 = t2[_n-1]
(1 missing value generated)

. replace t11 = 0 if[year < 1962]
(34 real changes made)

. replace t21 = 0 if[year < 1962]
(34 real changes made)

. gen y0lt = y0 * log(t1)
(75 missing values generated)

. gen y0lt1 = y0lt[_n-1]
(76 missing values generated)

. replace y0lt = y0 if[y0lt==.]
(24 real changes made)

. replace y0lt1 = y0 if[y0lt1 == .]
(30 real changes made)

. 
. nl (y = {d1}*aussie +{d2}*aust+{d3}*canada+{d5}*finland+{d6}*france+{d8}*iceland+{d9}*ireland+{d10}
> *japan+{d12}*norway+{d14}*spain+{d16}*uk+{d17}*us + ((1-{a2=.14}) + {a2=.14}*exp({a3=-2.3}*t1))*y1 
> + (1 - ((1-{a2=.14}) + {a2=.14}*exp({a3=-2.3}*t1)))*({b1=.4} + {b2=.7}*exp({b3=-.004}*t2)*y0) - {d=
> .13}*(y1-({b1=.4} + {b2=.7}*exp({b3=-.004}*t2)*y0 )) +{d=.13}*((1-{a2=.14}) + {a2=.14}*exp({a3=-2.3
> }*t11))*(y2-({b1=.4} + {b2=.7}*exp({b3=-.004}*t21)*y0))) if[year > 1960 & sample == 1], nolog vce(h
> c3)
(obs = 564)

Nonlinear regression                                 Number of obs =       564
                                                     R-squared     =    0.9993
                                                     Adj R-squared =    0.9993
                                                     Root MSE      =  .0196606
                                                     Res. dev.     = -2849.798

------------------------------------------------------------------------------
             |             Robust HC3
           y |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         /d1 |  -.0070707   .0062644    -1.13   0.260     -.019376    .0052347
         /d2 |  -.0157728   .0058984    -2.67   0.008    -.0273592   -.0041865
         /d3 |   .0127415   .0060223     2.12   0.035     .0009118    .0245711
         /d5 |   .0055211    .005236     1.05   0.292    -.0047641    .0158064
         /d6 |    .000426   .0049655     0.09   0.932    -.0093278    .0101799
         /d8 |   .0020014   .0068101     0.29   0.769    -.0113757    .0153786
         /d9 |  -.0171197   .0073113    -2.34   0.020    -.0314815   -.0027579
        /d10 |  -.0011263   .0063372    -0.18   0.859    -.0135746     .011322
        /d12 |  -.0035123   .0059311    -0.59   0.554    -.0151629    .0081383
        /d14 |  -.0035478   .0063935    -0.55   0.579    -.0161066    .0090111
        /d16 |  -.0150695   .0068655    -2.19   0.029    -.0285555   -.0015835
        /d17 |  -.0138418   .0094726    -1.46   0.145    -.0324489    .0047653
         /a2 |    .138316   .0299295     4.62   0.000      .079525    .1971071
         /a3 |  -2.097478   1.052936    -1.99   0.047     -4.16578   -.0291756
         /b1 |   .3802871   .0806819     4.71   0.000     .2218022    .5387721
         /b2 |   .7521068   .1360565     5.53   0.000     .4848485    1.019365
         /b3 |  -.0428131   .0150095    -2.85   0.005    -.0722965   -.0133297
          /d |   .1303669   .0756181     1.72   0.085     -.018171    .2789049
------------------------------------------------------------------------------

. display e(ll)
1424.8992

. 
. nl (y = {d1}*aussie +{d2}*aust+{d3}*canada+{d5}*finland+{d6}*france+{d8}*iceland+{d9}*ireland+{d10}
> *japan+{d12}*norway+{d14}*spain+{d16}*uk+{d17}*us + ((1-{a2=.16}) + {a2=.16}*exp({a3=-1.8}*t1))*y1 
> + (1 - ((1-{a2=.16}) + {a2=.16}*exp({a3=-1.8}*t1)))*({b1=.5} + ({b2=.54}+{b2us=-.95}*us)*exp({b3=-.
> 006}*t2)*y0) - {d=.13}*(y1-({b1=.5} + ({b2=.54}+{b2us=-.95}*us)*exp({b3=-.006}*t2)*y0 )) +{d=.13}*(
> (1-{a2=.16}) + {a2=.16}*exp({a3=-1.8}*t11))*(y2-({b1=.5} + ({b2=.54}+{b2us=-.95}*us)*exp({b3=-.006}
> *t21)*y0))) if[year > 1960 & sample == 1], nolog vce(hc3)
(obs = 564)

Nonlinear regression                                 Number of obs =       564
                                                     R-squared     =    0.9993
                                                     Adj R-squared =    0.9993
                                                     Root MSE      =  .0196277
                                                     Res. dev.     = -2852.723

------------------------------------------------------------------------------
             |             Robust HC3
           y |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         /d1 |  -.0126207   .0071263    -1.77   0.077    -.0266191    .0013778
         /d2 |   -.015544   .0057133    -2.72   0.007    -.0267667   -.0043212
         /d3 |   .0063829   .0065765     0.97   0.332    -.0065355    .0193012
         /d5 |   .0023935    .005479     0.44   0.662     -.008369    .0131561
         /d6 |  -.0003462   .0049228    -0.07   0.944    -.0100162    .0093238
         /d8 |   .0030647   .0067601     0.45   0.650    -.0102144    .0163438
         /d9 |    -.01492   .0068251    -2.19   0.029    -.0283267   -.0015133
        /d10 |  -.0026742   .0063926    -0.42   0.676    -.0152314    .0098831
        /d12 |   .0006502   .0053085     0.12   0.903    -.0097774    .0110777
        /d14 |  -.0056941   .0064435    -0.88   0.377    -.0183513     .006963
        /d16 |  -.0093694   .0060876    -1.54   0.124    -.0213274    .0025886
        /d17 |    .000863   .0088635     0.10   0.922    -.0165477    .0182737
         /a2 |   .1563556   .0390219     4.01   0.000     .0797039    .2330073
         /a3 |  -1.752866   .8693444    -2.02   0.044    -3.460541   -.0451895
         /b1 |   .5223765   .0847669     6.16   0.000     .3558668    .6888863
         /b2 |   .5305536     .13291     3.99   0.000     .2694751    .7916322
       /b2us |  -1.005902   .3813814    -2.64   0.009     -1.75506   -.2567448
         /b3 |  -.0625269   .0233731    -2.68   0.008    -.1084393   -.0166145
          /d |   .1271589    .076478     1.66   0.097    -.0230689    .2773868
------------------------------------------------------------------------------

. display e(ll)
1426.3613

. 
. nl (y = {d1}*aussie +{d2}*aust+{d3}*canada+{d5}*finland+{d6}*france+{d8}*iceland+{d9}*ireland+{d10}
> *japan+{d12}*norway+{d14}*spain+{d16}*uk+{d17}*us + ((1-{a2=.16}) + {a2=.16}*exp({a3=-1.8}*t1))*y1 
> + (1 - ((1-{a2=.16}) + {a2=1.8}*exp({a3=-.16}*t1)))*({b1=.5} + {b1us}*us + ({b2=.54}+{b2us=-.95}*us
> )*exp({b3=-.006}*t2)*y0) - {d=.13}*(y1-({b1=.5} + {b1us}*us + ({b2=.54}+{b2us=-.95}*us)*exp({b3=-.0
> 06}*t2)*y0 )) +{d=.13}*((1-{a2=.16}) + {a2=.16}*exp({a3=-1.8}*t11))*(y2-({b1=.5} + {b1us}*us + ({b2
> =.54}+{b2us=-.95}*us)*exp({b3=-.006}*t21)*y0))) if[year > 1960 & sample == 1], nolog vce(hc3)
(obs = 564)

Nonlinear regression                                 Number of obs =       564
                                                     R-squared     =    0.9993
                                                     Adj R-squared =    0.9993
                                                     Root MSE      =  .0196457
                                                     Res. dev.     = -2852.723

------------------------------------------------------------------------------
             |             Robust HC3
           y |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         /d1 |  -.0126756    .007456    -1.70   0.090    -.0273217    .0019705
         /d2 |  -.0155393   .0056921    -2.73   0.007    -.0267205   -.0043581
         /d3 |   .0062995   .0074266     0.85   0.397    -.0082888    .0208879
         /d5 |    .002347   .0058531     0.40   0.689    -.0091506    .0138445
         /d6 |  -.0003665    .005025    -0.07   0.942    -.0102373    .0095044
         /d8 |   .0030558   .0067686     0.45   0.652      -.01024    .0163516
         /d9 |  -.0148942   .0068826    -2.16   0.031     -.028414   -.0013744
        /d10 |  -.0027006   .0065494    -0.41   0.680    -.0155658    .0101646
        /d12 |   .0006774   .0053126     0.13   0.899    -.0097583    .0111131
        /d14 |  -.0057256   .0066438    -0.86   0.389    -.0187763    .0073251
        /d16 |  -.0093171   .0062504    -1.49   0.137     -.021595    .0029609
        /d17 |   .0008185   .0083933     0.10   0.922    -.0156687    .0173058
         /a2 |   .1562897   .0390456     4.00   0.000      .079591    .2329884
         /a3 |  -1.754448   .8700826    -2.02   0.044    -3.463581   -.0453151
         /b1 |   .5237816   .0997141     5.25   0.000     .3279098    .7196534
       /b1us |  -.0068895    .217362    -0.03   0.975    -.4338612    .4200821
         /b2 |   .5286102    .150869     3.50   0.000     .2322529    .8249674
       /b2us |  -.9750378   .8771351    -1.11   0.267    -2.698024    .7479487
         /b3 |  -.0628618   .0257812    -2.44   0.015    -.1135048   -.0122189
          /d |   .1271862   .0764717     1.66   0.097    -.0230298    .2774022
------------------------------------------------------------------------------

. display e(ll)
1426.3616

. log close
      name:  <unnamed>
       log:  /users/jjacksn/documents/pathdep/conference/pdtable1.log
  log type:  text
 closed on:  15 Jun 2011, 21:18:12
-----------------------------------------------------------------------------------------------------
